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Abstract 

Background: The absence of horns, called polled phenotype, is the favored trait in modern cattle husbandry. To 
date, polled cattle are obtained primarily by dehorning calves. Dehorning is a practice that raises animal welfare 
issues, which can be addressed by selecting for genetically hornless cattle. In the past 20 years, there have been 
many studies worldwide to identify unique genetic markers in complete association with the polled trait in cattle 
and recently, two different alleles at the POLLED locus, both resulting in the absence of horns, were reported: 
(1) the Celtic allele, which is responsible for the polled phenotype in most breeds and for which a single candidate 
mutation was detected and (2) the Friesian allele, which is responsible for the polled phenotype predominantly in 
the Holstein-Friesian breed and in a few other breeds, but for which five candidate mutations were identified in a 
260-kb haplotype. Further studies based on genome-wide sequencing and high-density SNP (single nucleotide 
polymorphism) genotyping confirmed the existence of the Celtic and Friesian variants and narrowed down the 
causal Friesian haplotype to an interval of 145 kb. 

Results: Almost 6000 animals were genetically tested for the polled trait and we detected a recombinant animal 
which enabled us to reduce the Friesian POLLED haplotype to a single causal mutation, namely a 80-kb duplication. 
Moreover, our results clearly disagree with the recently reported perfect co-segregation of the POLLED mutation 
and a SNP at position 1 390 292 bp on bovine chromosome 1 in the Holstein-Friesian population. 

Conclusion: We conclude that the 80-kb duplication, as the only remaining variant within the shortened Friesian 
haplotype, represents the most likely causal mutation for the polled phenotype of Friesian origin. 



Background 

The bovine polled phenotype, i.e. the absence of horns, has 
huge practical importance for breeders and is of special 
biological interest to geneticists. Currently, the world cattle 
population is estimated to be about 1.3 billion heads [1] of 
which a large proportion is horned. Although historical re- 
cords indicate the presence of naturally polled cattle as far 
back as ancient Egypt [2,3], until recently, horned cattle 
remained the desired phenotype because it simplified teth- 
ering and attachment to harnesses. However, in modern 
husbandry systems, such practices are no longer used and 
the presence of horns increases the risk of injuries to han- 
dlers and animals, especially as the housing densities have 
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increased. Thus, the presence of horns induces consider- 
able economic losses in the cattle industry due to dehorn- 
ing practices and the treatment of subsequent secondary 
infections, but also due to carcass and leather deterioration 
[2,4-7]. Today, dehorning of cattle at an early age is com- 
mon and although records are not exhaustive, there are 
probably more than 100 million calves that are dehorned 
each year. Since all dehorning methods are invasive and 
thus raise animal welfare issues, the possibility of breeding 
genetically polled cattle is a promising alternative [4,8]. In 
addition, understanding the molecular mechanisms that 
are involved in the inhibition of the bovine horn bud dif- 
ferentiation will contribute knowledge on the mechanisms 
that underlie ectopic expression in mammals [2], 

Since 1906, the polled phenotype is known to be inher- 
ited as an autosomal dominant trait [9] and the POLLED 
locus was mapped to bovine chromosome 1 (BTA1 for Bos 
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taurus) in 1993 [10]. In the past years, the position of 
the POLLED locus on BTA1 was refined [11,12] and can- 
didate causal mutations were identified [2,4,5]. Recently, 
Medugorac et al. [4] demonstrated the existence of at 
least two different alleles at the POLLED locus in cattle. 
While a complex 202-bp insertion-deletion (InDel), referred 
to as P202id> was identified in various cattle breeds from 
Scandinavia, Scotland, England, the Channel Islands and 
France down to the Alpine region (thus also referred to 
as the Celtic allele P c )> a 260-kb haplotype including five 
candidate mutations (three SNP: P G i654405A> Pci655463T> 
Pci76858?a; two InDel: P 5ID , Pgokbro) was reported to be per- 
fectly associated with the polled phenotype in the Holstein- 
Friesian breed (thus also referred to as the Friesian allele 
P F ) [4]. Further studies based on genome-wide sequencing 
and high-density SNP genotyping confirmed the existence 
of these Celtic and Friesian variants and narrowed the 
causal Friesian haplotype down to an interval of 145 kb [2] 
that included only two of the five initial candidate muta- 
tions i.e. SNP P C i768587A and InDel P 80 kbn> However, high- 
throughput sequencing of this 145-kb interval identified 
two new candidate variants i.e. SNPs Pti764239C and 
Pgi855898A> and genotyping of an enlarged cattle breed 
panel revealed the presence of the Celtic allele in new 
breeds from Iceland, Scandinavia, UK, Ireland and France 
[2]. In 2013, Allais-Bonnet et al. [2] provided evidence that 
introgression of the Friesian allele was ongoing in some 
French and German cattle breeds for which the Celtic al- 
lele was predominant and Glatzer et al. [5] reported that 
SNP AC000158:g.l390292G > A, hereinafter referred to as 
rsl34968598, co-segregated perfectly with the POLLED 
locus in a set of 443 Holstein-Friesian animals. This SNP 
is located within intron 3 of the interferon gamma receptor 
2 (IFNGR2) gene and is close (1113 bp) to a SNP for 
which the most significant selection signature was ob- 
served in a group of polled cattle [13]. However, these two 
SNPs are located 378 kb from the above mentioned 
Friesian haplotype. 

This study aimed at (1) refining the 260-kb haplotype 
which, based on previous results, is associated with the 
polled phenotype of Friesian origin [4] and (2) examining 
the co-segregation of SNP rsl34968598 and the POLLED 
locus in Holstein-Friesian cattle. These aims are highly in- 
teresting from a practical and economic point of view for 
the most important global cattle breed Holstein-Friesian. 
Identification of the causal candidate mutation in the 
Friesian haplotype will also contribute to decipher the 
molecular mechanisms that underlie this complex de- 
velopmental phenotype with allelic heterogeneity. 

Methods 

Ethical statement 

No formal ethical approval was required, since no new 
tissue samples were collected for this study. Indeed, for 



most animals, the DNA samples from previous studies 
[4,12] were in stock in our laboratory and the remaining 
DNA samples were obtained from routine testing (paternity 
and polled phenotype test) by appropriate cattle breeding 
organizations. 

Animal samples 

A total of 5993 bovine individuals were genotyped for the 
Celtic allele P2021D and for the two InDel mutations that 
flank the Friesian haplotype i.e. P 5ID and PgokbiD (Figure 1A) 
using previously described methods [4]. A small propor- 
tion of these animals was also genotyped for the three 
other mutations in the Friesian haplotype i.e. (i) animals 
that were shown to carry both the Celtic and the Friesian 
alleles on two different chromosomes (-1.2%), (ii) animals 
that were recombinant within the Friesian haplotype 
(-0.02%) and (iii) a few animals that displayed an initial 
ambiguous genotyping result (-1.5%). To refine the Frie- 
sian haplotype, we included the two recombinant animals 
and all their available ancestors and relatives, which 
constituted a pedigree of 13 genotyped Holstein-Friesian 
animals (Figure IB). 

To test the previously reported [5] co-segregation of 
the POLLED locus with SNP rsl34968598, which is lo- 
cated outside of the 260-kb haplotype, we used 130 ani- 
mals that were representative of the global polled and 
horned Holstein-Friesian cattle population. These 130 
animals included animals that we also used to refine the 
Friesian haplotype and bulls that had already been ana- 
lyzed in our previous study [4]. Of the 130 animals, 69 
were horned (pp), 54 were heterozygous polled (Pp) and 
seven were homozygous polled (PP). The POLLED geno- 
type of each polled bull based on progeny testing was 
available from the respective breeding organizations. 
The progeny test considered phenotypes of at least nine 
progeny per bull mated to horned dams. 

Sequencing of the Pji 764239c candidate mutation 

In our previous study, the region captured around the SNP 
Pti7642390 one of the two new candidate SNPs reported 
by [2], was sequenced by high- throughput sequencing, but 
the sequence coverage was low [4]. Therefore, this region 
was re-sequenced by conventional Sanger sequencing [14] 
using DNA samples from four horned animals for which, 
in the previous study, the sequence capture data suggested 
a discrepancy between the known phenotype and the ex- 
pected genotype, although the coverage was too low to 
infer reliable genotypes. 

Refining the Friesian haplotype 

In order to refine the 260-kb haplotype that perfectly asso- 
ciates with the POLLED locus of Friesian origin, we geno- 
typed 13 Holstein-Friesian individuals (that formed a small 
pedigree with two recombinant bulls and 11 relatives) for 
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Figure 1 Polled phenotype of Friesian origin. A) Relative position of the candidate mutations for the polled phenotype of Celtic (P c ) and 
Friesian (P F ) origin. All reference sequence (RefSeq) alleles are shown in blue and POLLED candidate mutations in orange. The insertion-deletion 
events (P 5 | D , P 20 2id and P 8 okbiD) are presented as bars and SNPs as lines. Note that this figure does not include the fairly distant SNP rsl 34968598 
at position 1 390 292 bp. B) Pedigree of 13 Holstein-Friesian animals used to narrow down the interval containing the POLLED mutation of 
Friesian origin. Color code is as above. The recombinant haplotypes of the ROSSI and ELBOY RED bulls are shown with the corresponding colors. 
Polled homozygous (PP) individuals are represented by solid circles (females) and squares (males); declared polled heterozygous animals (Pp) by 
half-filled symbols; horned animals (pp) by empty symbols; missing and ungenotyped individuals are marked with a diagonal line. 



the Celtic allele P2021D and the four Friesian variants P 5ID , 
Pgi654405A> Pci768587A and P 80 kbiD as described in 
Medugorac et al [4], Furthermore, of the two additional 
Friesian SNPs P G i855898A and P T i764239C reported in 
Allais-Bonnet et al [2], we included only the candidate 
causal mutation Pgi855898A because preliminary Sanger 
sequencing results indicated that Pti764239C segregated in 
horned animals. Thus, the 13 Holstein-Friesian individuals 
were also genotyped for Pgi855898a using PCR-RFLP. 
Therefore, PCR products of the region around Pgi855898A 
were amplified using the primers reported in [2], digested 



with the restriction enzyme PstI and the digestion 
products were size-separated by electrophoresis on a 
2% ethidium-bromide stained agarose gel. 

Genotyping of SNP rsl 34968598 to test co-segregation 
with the POLLED locus 

The 130 animals, including those used to refine the haplo- 
type, were genotyped for the recently published mutation 
within the gene IFNGR2 (rsl34968598) as described in 
Glatzer et al. [5]. 
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Results and discussion 

Genotyping of SNP rsl 34968598 to test co-segregation 
with the POLLED locus 

Recently, Glatzer et al [5] reported that SNP rsl34968598 
co-segregated with the POLLED locus in a set of 443 (15 
PP, 71 Pp, 357 pp) Holstein-Friesian animals [5]. In our 
study, genotyping of 130 animals representative of the glo- 
bal polled and horned Holstein-Friesian cattle population 
(7 PP, 54 Pp, 69 pp) for SNP rsl34968598 revealed dis- 
crepancies between the progeny-based genotypes and the 
SNP-based genotypes at the POLLED locus in seven cases 
(see Table 1). Of the 69 genotyped horned individuals 
(POLLED genotype pp), 67 were homozygous G/G for 
rsl 34968598 as expected according to [5] but two bulls 
were heterozygous A/G. Of the 54 proven heterozygous 
polled (Pp) bulls, 51 were heterozygous A/G and three 
were homozygous G/G. Finally, five of the seven proven 
homozygous polled bulls (PP) were homozygous A/A, but 
two were heterozygous A/G. In conclusion, -5% of the 
genotypes at SNP rsl34968598 are inconsistent with the 
previously reported perfect co-segregation between this 
SNP and the POLLED locus in Holstein-Friesian cattle. 
Similar to Glatzer et al. [5], our sample of 130 animals was 
representative of the global Holstein-Friesian population 
(i.e. from Canada, the Czech Republic, Denmark, France, 
Germany, Italy, the Netherlands and the United States of 
America) and among the polled individuals that exclude 
SNP rsl34968598, some are major breeding bulls like the 
homozygous polled PREMIUM PP bull which has been 
extensively used in recent years. 

Refining the Friesian haplotype 

After our genotyping data excluded the rsl34968598 
SNP on BTA1 at position 1 390 292 bp as a candidate 
mutation, we re-analysed our sequence capture data and 
partly re-sequenced a region around the Pxi764239C can- 
didate mutation of the 145-kb haplotype of Friesian ori- 
gin reported in Allais-Bonnet et al. [2]. Several horned 
animals were found heterozygous C/T at SNP Pxi764239C 
and thus this SNP was also excluded as a candidate mu- 
tation. In order to refine the Friesian haplotype, which 

Table 1 Seven bulls for which the progeny-based POLLED 
genotype and the genotype at SNP rsl 34968598 did not 



match 


Breeding bull 


International ID 


Polled 


rsl 34968598 


ELBOY RED 


FRAM0001 09031 700 


PP 


A/G 


ROSSI 


NLDM000538994750 


PP 


A/G 


JOSEF PS 


DEUM000768421954 


Pp 


G/G 


LEGOLAS P 


DEUM000352976294 


Pp 


G/G 


PRIMUS P 


DEUM000536884209 


Pp 


G/G 


PABLO PP 


DEUM000769554302 


PP 


A/G 


PREMIUM PP 


USAM0001 3941 9625 


PP 


A/G 



contains the causal mutation for the polled trait in most 
Holstein-Friesian breeding lines, a small pedigree includ- 
ing two recombinant bulls and relatives (Figure IB) was 
genotyped for the six known polymorphisms in this re- 
gion P202ID, P5ID> PG1654405A> Pci768587A> Pa1855898G and 

P 8 okbiD (for physical positions (according to UMD3.1) 
see Figure 1A). 

ELBOY RED, a horned male offspring of LAWN BOY P- 
RED, is a recombinant between P G i654405A and P C i768587A, 
which excludes P 5ID and Pgi654405A as well as the fairly dis- 
tant rsl34968598 SNP as candidate mutations for the 
polled phenotype. Previously, we reported and confirmed 
this recombination in 20 horned offspring of ELBOY RED 
[2]. A second horned offspring of LAWN BOY P-RED, the 
recombinant bull ROSSI, carries a large region of the Frie- 
sian haplotype that is associated with the polled trait and 

that Comprises P 5ID , P G i 6 54405A> Pci768587A and P G 1855898A 

(Figure 1A and B). The genotyping data of ROSSI located 
the recombination event in an interval of 53 kb between 
Pgi855898A and PgokbiD (beginning at g. 1909352) and clearly 
excluded P 5ID , P G i654405A> Pci768587A and P G i855898A as 
candidates for the polled phenotype of Friesian origin. 

Finally, based on the combined sequencing and genotyp- 
ing data from our previous studies [2,4] and this study, 
the only remaining causal candidate within the Friesian 
POLLED haplotype is the 80-kb duplication PgokbiD- This 
conclusion is inconsistent with that of Glatzer et al. [5] 
who reported the existence of a progeny-proven heterozy- 
gous polled Holstein-Friesian sire that carried neither the 
PsokbiD mutation nor a Celtic mutation. Unfortunately, 
since the identity of this bull was not provided, we could 
not verify these results. However, our personal experience 
from different previous analyses that we carried out is, 
that reliable genotyping of P2021D and PgokbiD is not trivial 
and that the genotyping-procedure should be validated 
using anonymous samples. Thus, there is a possibility that 
the bull reported in [5] does carry either P2021D or Pgoi<bn> 
In addition, even if we had demonstrated the absence of 
both P2021D and PgokbiD i n the genome of this buU further 
analyses would still have been necessary to show that it 
does not carry another spontaneous mutation, which 
could be responsible for a dominant horn defect syn- 
drome. Such mutations have been reported to occur and 
can be confused with the regular POLLED locus [2,15,16]. 
Therefore, the existence of a new candidate mutation at 
the POLLED locus of this animal should first be investi- 
gated and supported by linkage analyses on its relatives 
before refuting the highly probable causality of PgokbiD for 
the polled phenotype in Holstein-Friesian cattle. 

Conclusions 

Based on the genotypes of almost 6000 animals, we identi- 
fied two bulls with a recombination within the 260-kb 
Friesian haplotype and narrowed this haplotype down to 
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an interval of 53 kb on BTAL Moreover, the results from 
these two recombinant bulls and related individuals dis- 
agree with the perfect association between the POLLED 
locus and a SNP within intron 3 of the IFNGR2 gene previ- 
ously reported by Glatzer et al [5]. Combining the ex- 
haustive sequencing data and high-density SNP genotyping 
results from previous studies [2,4] and those reported here, 
we provide evidence that supports the PgokbiD mutation as 
the only remaining and most likely causal mutation for the 
polled phenotype of Friesian origin. In addition to Pgokbn> 
P2021D (the POLLED mutation of Celtic origin) segregates 
in some Holstein-Friesian families and sporadically com- 
bines in some homozygous polled but heterogeneous (P c / 
P F ) animals (Figure IB). Neither PgokbiD nor P2021D are 
located in any known coding regions [2,4] but these muta- 
tions, even as single copy, result in the inhibition of horn 
development at a very early stage in cattle [2,4] . Therefore, 
future investigations of the molecular mechanisms that 
underlie the inhibition of horn development in cattle will 
enrich our knowledge not only on the mechanisms that 
are involved in bovine horn bud differentiation but, also, 
more generally, on ectopic expression in mammals. 
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